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Abstract 

Using the Theory of Scattering in Restricted Geometries developed by A. Lupu-Sax as a 
starting point, we present a comprehensive multi-channel theory of atom-atom scattering 
in tight atom waveguides. 
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Conventions and notations 

Square root conventions 

Throughout the text we will be using two complementary conventions for the complex square root 
function: the usual convention \fz defined as 

^J\z\e i( t> = y/\z\e i<i>/3 < (f> < 2tt , 
and a non-traditional \fz one defined as 

= y/\7\e^ /2 ~ 2tv < < . 



Some notations 

Gjj(E) Green's function of a Hamiltonian H 

Tfi v(E) T-matrix of a perturbation V to a Hamiltonian H 

G(E) Green's function of the harmonic waveguide with no scatterer present 

Xjj(E) regular part of the Green's function of a Hamiltonian H at the origin 

x{E) regular part of the harmonic waveguide Green's function at the origin 

£ = E/2Tilu±_ — 1/2 rescaled and renormalized energy 

a±_ = y/h/(j,ui± transverse harmonic oscillator length 

\nm) eigenstates of the two-dimensional harmonic oscillator 

a three-dimensional scattering length 

g — 2Tvh 2 a/fi three-dimensional coupling constant 

OiD one- dimensional scattering length 

gin = —h 2 /[J.aiD one- dimensional coupling constant 



1 INTRODUCTION 



Atom waveguides are a fundamental component of Atom Optics, and are expected to play 
an important role in Atom Interferometry and Quantum Computing applications. To ensure 
proper coherence as atomic beams propagate through the waveguides, effort should be made to 
avoid decoherence-inducing mechanisms such as collisional losses and collisional phase shifts. 
This clearly requires a detailed understanding of the effects of quasi-one-dimensional confine- 
ment on atom-atom collisions. In particular, we will see that in the few-mode regime, necessary 
for coherent propagation, the free-space estimates of the collisional effects are no longer valid 
and a waveguide-specific theory is needed. 

Such one-dimensional interacting atomic quantum gases have recently been attracting sig- 
nificant theoretical and experimental interest, having become accessible via adiabatic transfer 
from atomic Bose condensates to highly elongated tight cigar-shape traps. Here the question of 
effective one- dimensional coupling constants becomes important, as they play a crucial role in 
determining whether or not the ground state posses long-range order (coherence). In the regime 
of tight confinement it is the virtual excitation of the transverse modes which play a crucial 
role, leading to the confinement-induced renormalization of interatomic interactions 0E] (see 
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also the two-dimensional analog of such a renormalization described in [3] and the fcrmionic 
p-wave analog in pi]). 

In this paper we present a comprehensive theory of atom-atom scattering in atom waveg- 
uides. Our theory correctly takes into account both quantization of the transverse motion 
and the transverse renormalization of the collisional strength, otherwise inaccessible by the 
free-space scattering theory. 

From the formal point of view the problem reduces to a Schrodinger equation for two 
atoms in a harmonic waveguide. In this paper we concentrate on the universal properties of 
the waveguide scattering which are governed uniquely by the scattering length. The Theory of 
Scattering in Restricted Geometries developed by A. Lupu-Sax allows us to describe these 
properties without invoking the full interaction potential, which is very often unknown. 



2 FORMULATION OF THE SCATTERING PROBLEM 



We begin from the Hamiltonian for two atoms under transverse harmonic confinement and 
subject to an arbitrary interaction potential 

tt ^ t-,2 ft t-,2 ,1 2 2,1 22 

+ V(ri-r a ) (2.1) 

where m% and ra-i are the atomic masses, u± is the transverse trap frequency, and V 2 and 
Tix are the Laplacian and radial coordinate of the i th atom, respectively. This Hamiltonian 
is separable in relative and center-of-mass coordinates R = (miri 4- m-2X , i)lM, r = ri — r2, 
M = mi + 77i2 being the total mass, yielding if 2 = H re i + Hcom, where 

h 2 1 

Hrei = -^V r 2 + -imjItI + V(r), (2.2) 

and 2 

Hcom = ~^%. + \MJiB? ± , (2.3) 

where \i = m\m-zl {m\ + 7712) is the reduced mass, and rj_ and Rj_ are the relative and center- 
of-mass radial coordinates, respectively. The center-of-mass Hamiltonian is that of a simple 
harmonic oscillator whose solution is known, hence we focus only on the relative motion of the 
two particles. This reduces the problem to a single particle of mass /x, subject to transverse 
harmonic confinement, which is scattered by an external potential V(r). The central equation 
which must be solved is therefore Schrodinger 's equation for the state of relative motion of the 
two atoms 

E — H — V] \ip(E)} = (2.4) 

where if = H re i and V is the interatomic potential. Determining the eigenstates of this 
Hamiltonian, in particular within the s-wave scattering approximation for V(r), is the central 
goal of this paper. 



3 GREEN'S FUNCTION AND T-MATRIX FORMALISM 



3.1 Definitions and theorems 

In this section briefly review the T-matrix formulation of scattering theory, which provides a 
convenient framework for approaching the present problem. Let us first introduce the retarded 
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Green's function for a system with the Hamiltonian H and energy E 



Gjj(E) = lim {E + ie- H) 



(3.1) 



We then define the T-matrix at energy E of the scatter V in the presence of the background 
Hamiltonian H (Fig^,,b) in the usual manner as 



f„y(E) = [l-VG^E) 

OO 



V 



n=0 



v. 



(3.2) 



the summation form being valid provided that there are no difficulties with convergence. 
Two relations on which we will rely heavily are the Lippman-Schwinger relation 

G 6+V (E) = G A {E) + G A (E)f A #{E)G A (E), 



(3.3) 



which relates the full Green's function of the system H+V to the unperturbed Green's function 
Gfi(E) and the T-matrix, and the Lupu-Sax formula (FigQJ;) 



T k y{E)= 1-7WE) G A (E)-G A ,(E) T^E), 



(3.4) 



which relates the T-matrix of the scatter V in the background Hamiltonian H (FigEJa) to 
the T-matrix for the same scatter but in a different background Hamiltonian H' (FigQJi). 
Derivations for these expressions are given in Appendices IB . 21 and IB . 31 respectively. 



3.2 Scattering theory 

In the continuous part of the spectrum of the total Hamiltonian H + V its eigenstates can be 
expressed as a sum of an incident and a scattered wave according to 



mE)) = \ME)) + \ME)), 
where \tpo(E)}, the 'incident' state vector, satisfies 

G A 1 (E)\ME))=0, 
we can then express the Schrodinger equation for the total system as 



G n 1 {E)~V (\ME)) + \ME))) = 0. 



(3.5) 



(3.6) 



(3.7) 



This equation is readily solved for the scattered wave in terms of the unperturbed Green's 
function and the T-matrix (see Appendix IB. lfl . yielding 

\ME)) - G 6 (E)f 6 ^(E)\ME)), (3.8) 
which will serve as the basis for our treatment of the present scattering problem. 
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Fig. 1. An artist view on the Lupu-Sax theorem. Here G /ree = G„ (E), G bound - = G 6o+trbound . (E), 
T free = f^^E), T bound - = f 6o+(rboand . ^ (E) . In our problem H corresponds to the kinetic 
energy, (j bound - i s the waveguide potential, and V s is the interatomic interaction potential. 
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3.3 Bound T-matrix as a generator for bound state energies 

It follows from Eq. IjH.lJI that G^-(E) is diagonal in energy representation and therefore takes 
the form 



GhW-I^^Te^+I^J dE E -E' + i e' (3 - 9) 

n j L >3 

where the first term on the r.h.s. sums over the discrete portion of the spectrum of H and the 
second term sums over all continuous bands of the spectrum. The limit e — > + is implied. 
From this expression we see that in the discrete (bound-state) part of the spectrum the poles 
of the Green's function correspond to the bound state energies. 

Consider now the particular case where the background Green's function G^(E) has only 
a continuous spectrum bounded from below by E = 0, i.e. 

Now consider a system consisting of this background Hamiltonian H plus a scattering potential 
V, described by the T-matrix Tg- y(E). Inserting equation (|3.1U|) into the Lippman-Schwinger 
Equation (|3.3|) for the total Green's function gives 

n (E)- fV lE ' ){E ' 1 I r HdE'dE" \ E ')^'\THyiE)\E")(E"\ 
Gh + v{E) -Jo ^ E-E> + ie + J J dEdE (E-E> + ie)(E-E»+ie) (3 ' 11} 

If a bound state of the combined system H + V having energy E n < exists, then the full 
Green's function must take the form 

Since the background Green's function has no pole at E n , it follows that the T-matrix itself 
must have a singularity at E = E n . Finding the poles of the T-matrix below the lower bound 
of the spectrum of the background Hamiltonian therefore gives a method for determining the 
bound state energies of the system H + V. 



4 S-WAVE SCATTERING REGIME: THE REFERENCE T-MATRIX APPROACH 

At first glance it may seem that finding the T-matrix is no easier than a direct solving of 
the Schrodinger equation 1)2 .4|) . We will demonstrate, however, that the T-matrix formulation 
allows for a self-consistent description of the low-energy part of the spectrum that uses the 
free-space low-energy scattering properties of the interaction potential as the only input. In 
addition the low-energy (s-wave) limit is isolated to a single well-defined approximation without 
requiring the ad-hoc introduction of regularization via a pseudo-potential. In this section we 
first outline this self-consistent low-energy treatment. We then solve for the T-matrix using 
the standard Huang- Fermi pseudo-potential, showing that the pseudo-potential reproduces the 
exact result in this situation. 

Let the unperturbed Hamiltonian H be a Hamiltonian for a single nonrelativistic particle 
in presence of a trapping potential U: 



2V72 



rv 



2„ U ^ 



(■#) , (4.1) 
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Assume also that the particle is' perturbed' by a scatterer given by 

(r\V\iP)=V(r)(r\iP) (4.2) 

localized around r = 0. In what follows we will derive a low- energy approximation for the 
T-matrix of the scatterer V in presence of H. It is important to note that, by definition, the 
T-matrix acts only on eigenstates of the unperturbed Hamiltonian, which we can safely assume 
to be regular everywhere (this is of course a constraint on the properties of the unperturbed 
Hamiltonian). In this case the zero-range s-wave scattering limit does not require any regu- 
larization of the T-matrix. By making use of the Lupu-Sax formula l|3.4|l . we first derive the 
correct form of the T-matrix in the low-energy s-wave regime without the introduction of a 
regularized pseudo-potential. In the following section, however, we will see that the results we 
obtain are in agreement with the standard Huang-Fermi pseudopotential approach to s-wave 
scattering. 

We begin our derivation by first specifying a 'reference' background Hamiltonian H' as 



(r\H'\ 



2fi 



E 



<i#>. (4.3) 



This Hamiltonian is that of a free particle, but with an explicit energy dependence included so 
that the eigenstates have zero wavelength at all energies. We note that this reference Hamilto- 
nian agrees with the free-space Hamiltonian in the zero-energy limit. While this Hamiltonian 
may seem strange, it is a valid reference Hamiltonian which turns out to be useful because the 
resulting T-matrix is energy independent for any scattering potential. The Green's function 
for this Hamiltonian is given by 

as can be verified by direct substitution into [E — H'] G^,(E) = I. In turn the T-matrix of 
the interaction potential V in presence of H' is independent of energy and can therefore be 
expressed as 

(r\f A ,^(E)\r , }=gD(r,v'), (4.5) 
where the kernel D is defined as normalized to unity, 

J drdr'L>(r,r') = 1. (4.6) 

The normalization coefficient g is then related to the 3-dimcnsional scattering length a accord- 
ing to 

2nh 2 a 



(4.7) 



A* 

a relationship which is derived explicitly in Appendix IB. 41 

Imagine that the kernel D(r,r') is well localized within some radius R. In perturbative 
expansions at low energies this kernel only participates in convolutions with slow (as compare 
to R) functions, in which case it can be approximated by a (^-function, 

D(r,r') ps<y(r)£(r / ). (4.8) 
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This straightforward approximation is the key to the s-wave scattering approximation. As we 
demonstrate in detail in Appendix IB. 51 this effectively replaces the exact reference T- matrix 
by its long- wavelength limit, so that the reference T- matrix assumes the form 

(r\f A ,^E)\r') k ' k '^ 1/R 9 S(r)5(r'), (4.9) 

which is equivalent to 

T 6 ,#(E)=g\0)(Q\, (4.10) 

where |0) is the position eigenstate corresponding to the location of the scatterer. In expression 
H4.9JI k and k' refer to the wavevectors of any matrices which multiply the T-matrix from the 
left and right, respectively. 

If we now substitute the above expression for the reference T-matrix into the Lupu-Sax 
formula (|3.4|) for the T-matrix under the background Hamiltonian H we arrive at 

oo 

= [l- g{0\G & (E)\0) + giOlG^Emy 1 g\0){0\. 



Making use of Eq. I|4.4|l . we introduce the function x(£), defined as 



X (E) = Urn 

r^O 



(r\G A (E)\0} + —^ 



(4.11) 



(4.12) 



from which we obtain the following simple expression for the T-matrix of the scatterer V in 
presence of the trap: 

E<^h 2 /uR 2 a5(r) 
(r\T 69 (EM) « /M l-gxW m " ^ 

As we will explicitly demonstrate for the case of transverse harmonic confinement, the singu- 
larity in bound Green's function is the same as that in the free-space Green's function. Hence, 
x{E) is the value of the regular part of the bound Green's function at the origin. 

5 THE FUNCTION X (E) 

In this section we derive the bound Green's function of the waveguide, from which we obtain 
an analytic expression for regular part at the origin x(-E), which gives the energy dependence 
of the bound T-matrix. We begin by reviewing the eigenstates of the waveguide potential, 
which we then use as a basis for expanding the bound Green's function and obtaining \ (E). 

5.1 Eigenstates of the waveguide Hamiltonian 

The Hamiltonian for the relative motion of two atoms in a harmonic waveguide contains two 
parts, the longitudinal free Hamiltonian H z and the transverse confinement Hamiltonian H±, 

H = H Z + H X , (5.1) 

where 

<r|^> = "|^(rhA>, (5.2) 
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and 



(r|iTx|V) = 



if 



Id Id 2 



2p \dp 2 p dp p 



2 f)A2 



2 2 
-UUJ ± P 



(5.3) 



The eigenstates of the transverse Hamiltonian, denoted by \nm), are well known and satisfy 
both 

H±\n, m) = huij_(2n + \m\ + l)|ran), (5-4) 



as well as 



L z \nm) — hm\ 



(5.5) 



where L z is the operator for angular momentum along the z-axis. We note that in this rep- 
resentation the radial and azimuthal quantum numbers n and m independently assume the 
values 

7i = 0,1, 2,..., oo, (5.6) 



and 

m = G,±l,±2,...,±oo. 
The eigenfunctions are given in r-representation by 

-1/2 



(5.7) 



(pcf)\nm) 



Traj_(n + \m\)\ 



e-5(p/ Q -L) 2 ( /9 / a _ L )l™l e jm ^[l m l](p 2 /ai), (5. 



where 



aj_ = 



(5.9) 



is the transverse harmonic oscillator length. Lastly we note that the value of |n, m) at the 
origin is given by 

'n| nm ) = (5.10) 



which is independent of n. 



5.2 The Green's function for the relative motion of two particles in a harmonic waveguide 

The bound Green's function is the solution to the equation 



E — H z — H i 



G{E) = I, 



(5.11) 



where / is the identity matrix and for simplicity we have taken Gjj(E) <-* G(E). Expanding 
this equation onto the set of states 



\nmz) = \nm) <g> \z), 



(5.12) 



where |nm) is the eigenstate of H± given by (|5.8|l and \z) is the eigenstate of z, satisfying 
(z\z') = 5(z — z'), then gives 



h 2 d 2 

E-?uj x (2n+\m\ + l) + —-^ 



(nmz\G(E)\n'm' z') = {nmz\n'm'z') 

= S n ,n'S m , m 'S(z - z ). (5.13) 
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We proceed by making the Ansatz for the Green's function 

(nmz\G(E)\n'm'z') = S nyn ,S m . m >a nm e n " mlz ~ z ' 1 

= 6 n , n ,5 m , m ,a nm [e^^-^eiz - z') + e -^™(*-*') e ^' _ ^44) 

where Q(z) is the Heavyside step-function. Differentiating (|5.14|) twice with respect to z gives 

d 2 „ 

— ~{nmz\G{E)\n'rn! z 1 ) — —^ 2 m (nmz\G(E)\n'm'z') +2ij nm a nrn (nmz\n'm'z'), (5.15) 
oz z 

so the Eq. (|5.13|) becomes 



E- hu±(2n+ \m\ + 1) - 



2 - 2 1 ft 2 <v a 

nmz|G(£/)|n m 2; )-N (nmz|n m z ) — {nmz\n m z ) 



2(i. 



(5.16) 

This equation is satisfied provided that 



h 

and 



rim = ^lE-hLo^2n+\m\ + l)} (5.17) 



j 2 M 



a nm = -1-2 • (5.18) 



nm 



Equation (|5.17(l is quadratic, and therefore has in general two solutions. We can determine 
which solution is needed, however, from the causal nature of the retarded Green's function. 
We should choose the solutions which propagate outwards from z — z' . 
By introducing the dimensionless energy 

- E ' (5.19) 



2huj±_ 2 

and making use of the relation (|5.9|l . we find that retarded Green's function can be expressed 
as 

a e i 7[V^-¥ \ z ~ z '\ 

{nmz\G(£)\n'm'z') = -i—j , = & n ,n'8m,m'- (5.20) 

2h yf£- n -M 

With this expression for the bound Green's function we can now proceed to compute the bound 
T-matrix via equation (|4.13(l . 

5.3 The function x{£) for a point scatterer in presence of a wave-guide 

In order to evaluate x(£) we must first compute the matrix element (r|G(£)|0). Expanding 
this matrix element in terms of transverse eigenstates gives 

(r|G(£)|0) = (p<f>\nm)(nmz\G(£)\n'm'0)(n'm'0\Q), (5.21) 

nmn'm' 

which, with the help of Eqs. I|5. 11)11 and i|5.2L)|) . gives 

00 __2- kfrT^Z \z\ 

g(v\G(£)\0) = -v^a^(p^0) e 7/ c > ( 5 - 22 ) 
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where the modified square root ^/ is defined in the "Conventions and notations" table above. 
Now by inserting Eqs. (|5.22l) and (|4.7I) into Eq. 14.121) we arrive at the expression 



9X(£) 



lim 



,]T(W>|nO> 



/n-£\z\ 



(5.23) 



While both terms in Eq. I|5.23|l diverge in the limit r — > oo, we now demonstrate that their 
difference remains finite and leads to an expression for x(£) in terms of a generalized Zeta 
function. 

First we assume, without a proof, that the multi- variable limit in Eq. I|5.23|) does exist, i.e. 
the directional single-variable limits r = sn, s — ► exist for all directions n and they are all 
equal to each other. This assumption allows us to deal with the limit along the z-axis only: 



9X{£) = Ijm 



E : 

n=0 



/n—£\z\ 



(5.24) 



where we have used the identity (|5.1U|I . We proceed by first replacing the a/\r\ term in Eq. 
<|5.24[1 with an integral expression via the identity 



a 
M 



a 



dn ■ 



\Jx — E 



(5.25) 



This allows us to write 



q\{£) = lim lim 

a± |*|-.QJV-)-oo 



E 

n=0 



N 



/n=£|r| 



dn ■ 



(5.26) 



It is now tempting to interchange the limit signs and thus get rid of the coordinate dependence. 
In order to be able to do that one have to prove the uniformity with respect to \z\ of the TV — * oo 
convergence of the expression in the square brackets H(JV, \z\), i.e. to prove that for every e 
there exists N*, the same for all \z\, such that 5(JV, |z|)— limAr^oo 3(JV, \z\) < e for all N > N*. 
Such a proof does exist: provided its length we do not exhibit it here. 
We arrive at the following expression for %(£ ): 



9X(£) = J im 



N 1 



2yfN-£ 



(5.27) 



One can now make use of the following theorem involving the Hurwitz Zeta function, an 
analytic generalized Zeta function described in the mathematical literature [Hj, 



In particular, 



C(s, a) = lim 

N— >oo 



1 



t N 

E c ?) ^ n 

^n=0 v ' 



1 



1 



1 - s (N + a) s 



Re(s) > 0, -2tt < arg(n + a) < 0. 



(5.28) 



C(l/2,a)= lim 

N— >oo 



N i 

y^ 



- 2 VN + ■ 



(5.29) 
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While this expression, valid for any N, may be taken as a definition of the Hurwitz Zeta 
function, it does not constitute an efficient method for computation. Most symbolic math 
software packages will have efficient algorithms, however in using them care must be taken 
with regards to the branch cut for fractional s. With this definition and taking s — 1/2 we 
arrive at 

gx {E) = -— C(l/2,-£). (5.30) 

aj_ 

By substituting this expression into Eq. I|4.13(l we arrive at the final expression for the 
long-wavelength T-matrix in the waveguide: 



We note that in the definition Ij5.28|l of the Hurwitz Zeta function, there is an ambiguity in 
the sign of the square root. 

6 RESULTS 

6.1 Multi-channel scattering amplitudes and transition rates 

From Equations (|3.8[) and (|5.31|) we find that the scattered wavefunction takes the form 

(nmz\M£)) = Hfi T (0\M£))- (6-1) 

i + ^C(i/2,-£) 

Now the matrix element (nmz\G(£)\0) can be determined by by making use of Eqs. I|5.20|l and 
(|5.10() . yielding 

(nmz\G(£)\0) = f dz ' (nmz\G(S)\n'm'z'){n'm'z'\0) 

ri'rn' 

i-^— \f£—n\z\ 



A* c e 



-i n ^ 2 S m ,o 7= (6.2) 



Inserting this expression into Eq. H6.1|) then gives 

l + ^C(lA-f) 



(nmz\ip s (8)) = -i- v — T 5 m , e _ — (0^o(g))- (6.3) 



Let us now assume that the incident wave has the longitudinal wave vector k and the transverse 
quantum numbers n and m, according to 

(r\M£)) = ( P 4>\nmy k \ (6.4) 

where the relation 

/ a ± k\ 2 \m\ 

£= \—) +U+ — (6 ' 5) 

gives the dependence of the scaled energy £ on the incident wavevector k. This incident 
wave is nonzero at the origin only for m = 0, hence only incident waves with zero angular 
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momentum will scatter. The value of the m = incident wave at the origin conveniently takes 
the n-independent value, 

Wq(£)> = -fe 9 -. (6-6) 



/na± 



Assuming henceforth m = 0, we can now express Eq. (|6.3() as 



(n'm'z\M£)) 



<W, 



[^+c(i/2,-f)] vr 



(6.7) 



From this expression it follows that the full wavefunction of the relative motion takes the form 



(!#(£)) - £ (/^|n'0) k',„e lfcz + /(fc„, «- fc„)„^„e ; 



fe„/|z| 



(6.8) 



n'=0 



Here we have introduced the even-wave transversely inelastic scattering amplitudes f(k n 



given by 



f(kr, 



2/ 



"n in < — n 



a±k n 



^ + C(l/2,-(s^) -n) 



(6.9) 



and the outgoing wave vector for the mode |n'0) 



fe*,/ — 



/ a ± k r< 



o_l V V 2 



+ n — n 1 



(6.10) 



from which the desired scattering probabilities can be computed. 

One can now easily compute the elastic and inelastic transition probabilities for collisions 
under transverse harmonic confinement. We begin by considering the asymptotic forms of the 
total wavefunction given by Eq. (|6.8() . which are given by 



lim (ri0z\i(>(£)) 



S n >, n e lkz + 6 [£ - n'} f(k v 



v n in' * — n c 



lim {n0z\^(£)) = 6 [£ - n'] f(k n ' <— k n ) n <^ n e 



7.—^ — OC 



(6.11) 
(6.12) 



Because of energy conservation, an inelastic collision results in a change in the longitudinal 
momentum k — > k n , so that the introduction of inelastic transmission and reflection coeffi- 
cients must be based on conservation of total incident and outgoing probability current. For 
longitudinal plane waves the probability current is given by the amplitude squared times the 
velocity, which leads to the inelastic transmission and reflection coefficients 



R(k n 



v n Jn < — n 



v n Jn * — n 



Q[£- n'] 
Q[£- n'} 



l£- 


n> 




n 


l£- 


ri 



\$n' ,n ~f~ fifen' fenjn'*— n\ 

\f{k n f < kn 



^njn'i — n\ i 



(6.13) 
(6.14) 



which are readily evaluated with the help of Eq. (|6.9(l . The sum of the transmission and 
reflection coefficients of the particular channel gives the corresponding transition probability 
governing the population exchange between the transverse vibrational levels. 
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Substituting the expression l|6.9|) to the transmission and reflection coefficients above we 
obtain the following set of kinetic coefficients: 



W^n(k) =1- W^ n (k) 

2 7?(£)-l/V^ 



f£ t [2 ± + C (l/2,l-^)] 2 + ^(f) 



W n ,^ n {k) = Q[£ -ri] -Jf 1 



T t V£~~^ [f ± + C(l/2, 1 - 5£)Y + rf{£) 
T^ n (k) = W^ n (k)-R^ n (k) 



£ t [f ± + C(l/2,I-^)] 2 + ? ? 2 (f) 
(fc) = WW n (fc)/2 



(6.15) 
(6.16) 

(6.17) 

(6.18) 
(6.19) 
(6.20) 



where W^, n (k) is the probability that after a collision of two particles with the relative mo- 
mentum k and relative transverse excitation n the particles will remain in the same transverse 
state, this probability is a sum of transmission T<_ >n (fc) and reflection R^ n (k) probabilities, 
W^ n (k) is the total probability of changing the transverse state, W n '«_ n (fc) is the probabil- 
ity of transition to a particular transverse channel, the transmission T n ' and reflection 
Rn'*—n{k) probabilities in channel-changing collisions are equal to each other. Here 



£ = n + £ t 
£ t = (fc a± /2) 2 
6£ = £-[£} < 8£ < 1 
[£] 



l 



n'=0 



V£^ 



(6.21) 
(6.22) 
(6.23) 

(6.24) 



where [. . .] is the integral part sign. 

6.2 Single-channel scattering and effective one-dimensional interaction potential 

Let us now consider the special case of a single-channel scattering: 

< £ < 1, n = 

In this case we have 



a± aj_ 



9X(£) 

Using an alternative representation for C(l/2, 1 — £) 



«iai-*> + ^ 



«l/2,l-£)= Jim 



N l 
£ i/zr. 



V^£ 



2VN 



and making use of the expansion in powers of £ , 



1 



. AT ^ — ' 



(2j-l)" 5 ,- 



|£| < 1, n>0, 



(6.25) 



(6.26) 



(6.27) 



(6.28) 
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allows us to write 

C(l/2,l-f) = C(l/2) +£{£), (6.29) 

where 

3=1 J ' 

which clearly separates the zero energy limit from the finite energy corrections. 
According to (|6.9(l this leads to 

fe(k) = - T 1 - (6-31) 

where / e (fc) = f(ko <— fco)o^o is the even single-channel scattering amplitude and 



a 1D = - T 



— +C(l/2) (6.32) 



a 



with C(l/2) = —1.4603 . . ., is the effective one-dimensional scattering length which agrees with 
the result in pQ. 

It is now tempting to introduce an effective one-dimensional interaction potential in such a 
way that its scattering amplitude, introduced through the one-dimensional scattering solution 
as 

ip(z) z ~^°° exp(ikz) + f e (k) exp(ik\z\) + / (fc)sign(z) exp(ik\z\) , (6.33) 

matches (|6.31|l . i.e. solve the corresponding one-dimensional inverse scattering problem. It 
turns out that this problem is ill-posed due to the presence of open transverse channels unac- 
cessible within the one-dimensional model. Nevertheless one may pose the following problem: 
find a one-dimensional potential, whose scattering amplitude reproduces the exact one (|6.31|) 
with the relative error 0(k 3 ). Such an object does exist, and it is represented by a zero-range 
scatterer 

v(z) = g 1D S(z) (6.34) 

of a coupling strength gw = —h 2 /fiaio- 

Notice now the resonant behavior of g\o showing a Confinement Induced Resonance (CIR) 
at a — aj_/|C(l/2)|. The effect was recently interpreted in terms of Feshbach resonance between 
ground and excited vibrational manifolds. The resonance has been confirmed by numerical 
calculations with finite-range potentials, at both two-body |2] and many-body [7] level. 



6.3 Bound states 

As it has been discussed in the section 13.31 the poles of the full T- matrix of the problem (the 
T-matrix 15.311 in our case) correspond to the bound states. We get the following eigenvalue 
equation 

C(l/2, -£) = -— (6.35) 
a 

The detailed analysis, interpretation, and testing of the solutions against finite-range models 
is presented in [21 - It turns out that for any set of parameters there exists one and only one 
bound state. For small positive three-dimensional scattering length it converges to the free- 
space three-dimensional bound state. For small negative scattering length the bound state 
corresponds to the bound state of the one-dimensional potential (|6.34|l . 
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7 CONCLUDING REMARKS, RELATED WORKS, AND OPEN PROBLEMS 

In this paper we have demonstrated that the low-energy free-space properties of a scatterer 
are sufficient to describe its low-energy behavior in a non-free environment. The T- matrix for- 
malism and the Lupu-Sax connection formula in particular serve as a powerful bridge between 
the two. 

Several related works ought to be mentioned. Scattering of bosons in two-dimensional 
(planar) harmonic waveguides has been successfully treated in jjj- Fermions in a linear guide 
were considered in 0j via the K-matrix formalism and the corresponding one-dimensional 
amplitudes were explicitly computed. In an TV-body setting N — 2 particles can be interpreted 
as the background potential for a given pair [Hj. 

One would expect that the inverse scattering problem posed in section 16.21 can be solved 
with an accuracy higher than existing 0(k 3 ). The solution will allow to improve the accuracy of 
the many-body Monte-Carlo numerical models without introducing the transverse dimensions, 
otherwise making the computation harder. 

Shallow atomic guides with only a few transverse bound levels constitute a significant 
challenge. Unlike for harmonic guides the separation of relative and transverse degrees of 
freedom will be lifted leading to new collision channels. Presence of continuum spectrum for 
virtual transverse excitations may significantly enhance the renormalization effects or even lead 
to new resonances. 

APPENDIX 
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B APPENDICES 

B.l The T-matrix, the potential and the background Green's function 

In this section we derive various useful relationships between the Green's function of a given 
background Hamiltonian, G^(E), the scattering potential, V, and the T-matrix, T^y(E), 
which connects the incident wavevector \ipo(E)) with the scattered wave \4> S (E)) according to 
Eq. H3.8[) . By making use of the definition of the T-matrix (|3.8|l and the fact that the incident 
wave satisfies the unperturbed Hamiltonian (|3.6I) . we can cast Schrodinger's equation (|3.7[l in 
the following form 



t h,v( E ^- VG 6( e ) t h,v( E )- V \ME)) 



(B.l) 



Since this equation must be satisfied for any \ipo(E)) it follows that 



T A #(E) - VG 6 (E)T A ir (E) -V = 0. 



(B.2) 
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This equation can be readily solved for either the T-matrix V (E) or the potential V. Solving 
first for TW y(E) then gives 



TW£) = 1-VG A (E) 



V. 



(B.3) 



which gives an expression for the T-matrix at energy E in terms of the perturbation V and 
the background Green's function Gjj(E). We can also solve Eq. (|B.2|) for V, which yields 



V = f Av (E) l + G^ffjyiE) 



(B.4) 



or by starting from the Hermite conjugate of Eq. I|3.7(l we can similarly arrive at the equivalent 
expression 



V = 



l+f ^{E)G A {E) ffr v (E) 



(B.5) 



B.2 The Lippman-Schwinger equation 



In this section we derive the relation between the full Green's function of the Hamiltonian 
H + V in terms of the Green's function for the Hamiltonian H and the T-matrix. The full 
Green's function can be defined by 



which, after substituting Eq. (|B.4|) yields 



GZ" .(E) = G^ f Av (E) 1 + G 6 (E)f AtV (E) 



(B.6) 



(B.7) 



Multiplying from the left by G^ +V (E) and then from the right by 1 + Gg(E)ffy ^(E) G^(E) 
then yields the Lippman-Schwinger relation 

G 6+V (E) = G 6 (E) + G 6 (E)f 6v (E)G 6 (E) (B.8) 

which relates the full Green's function to the 'background' Green's function and T-matrix. 

B.3 The Lupu-Sax formula 

In many situations one may not have knowledge of the potential V, but rather have direct 
knowledge only of the T-matrix, Tfj, V (E) with respect to some background Hamiltonian H' . 
If additional external fields are applied, then ideally one would like an expression for the T- 
matrix of the same perturbation in the presence the new background Hamiltonian y (E) . 

We can derive the required expression directly from the relations l|B.4jl and (|B.5(I . Since 
V is equal to itself we can equate expression <|B.4() for V in terms of H' to the equivalent 
expression <|B.5(1 in terms of H, which yields the relation 



(B.9) 



l + t 6l p{E)& Al {E) f ft ,f(E)=ff tv (E) l + G^Ejf^yiE) 



Multiplying from the left by 1 + f 6 , V (E)G 6 ,(E) and from the right by 1 + G ij (E)f fl V (E) 
then gives 



f k y(E) l + G 6 (E)fz v (E) = l + f^ v (E)G 6 ,(E) f 6 *(E). 



(B.10) 
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We note that in this equation the unknown operator y (E) appears linearly and always on 
the far-right side. Hence, solving for V (E) is straightforward and results in the Lupu-Sax 
formula [H] relating the T-matrices of the same perturbation but in two different background 
Hamiltonians H and H'\ 

f &v {E)=[l-f &l y{E)(G ti {E)~G &l {E))Y\ k , v {E). (B.ll) 
B.4 The free-space scattering length 

In this section we will related the 3-dimensional scattering length to the normalization of the 
kernel of the reference T- matrix 1)4. 5fl . In Section 0] we introduced the reference Hamiltonian 
H' = El f ree + E, where Hf re e is the free-space Hamiltonian containing only a 3-dimensional 
kinetic energy term. We note that for E = the reference Hamiltonian and the free-space 
Hamiltonian agree, hence we can use the Green's function and T-matrix of H' to solve the 
free-space scattering problem at zero energy. The free-space solution for scattering from the 
potential V is therefore given at E = by 

\4> a (p))=G A , (0)4,^(0)1^0(0)) (B.12) 

where Gjj, (0) is given by Eq. 14.4f> and T^, y(0) by Eq. (|4.5|) . Expanding this expression onto 
position eigenstates and making use of (|4.4H and (|4.5I) gives 

(r|^(0)) = J dv'dv' ^^(r> (0)). (B.13) 

For zero energy the incident wave is given by (r|^/;o(0)) = 1, hence by making the assumption 
that D(v' , r") is localized at r', r" sa we find that the limit as r — » oo is given by 

lim(i# s (0)) = -J* 1 f dv'dv" D{v',v") 
r ^°° 2-kTi \v\ J 

= -^lA- (B-14) 
2-itTt 2 \v\ V 7 

Because the potential is assumed to vanish as r — > oo, we know that the asymptotic form 
of the scattered wavefunction is given by a spherical plane wave with zero kinetic energy, 
characterized by the zero-energy scattering amplitude / s (0), 

lim<r|^ s (0))=/ s (0)i (B.15) 



Comparison with Eq. (|B.14|l then shows that the zero-energy scattering amplitude is 

f s (E = 0)=-J^. (B.16) 
The asymptotic form of the total wavefunction is then 

lim(r|^(0)) = l--^i. (B.17) 



By definition, the scattering length is the radius of the first node of wavefunction at zero energy. 

(B.18) 



Setting Eq. I|B.17|) equal to zero and solving for r then gives 



2ttH z 

which relates the scattering length a and the normalization of the reference T-matrix, giving 
the standard expression g = 2nTi 2 a/ fj,. 
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B.5 The s-wave scattering approximation 

In this section we will discuss the long-wavelength properties of the reference T-matrix ((4.5(1 
and motivate the s-wave scattering approximation in which we replace this T-matrix by its 
long- wavelength (low energy) limit. From Eq. 1(4.5(1 we have 



(T\f 6 ,^\r')=gD(r,r'), 



(B.19) 



where the kernel D(r,r') is defined as normalized to unity. We note that the exact expression 
for D(r, r') via Eqs. 1(3.2(1 and 1(4.4(1 is given by 



gD(r,r') 



V(r)S(r - r 

3 



V(r)V(r' 



-/<■ 



2irh 2 



V{r)V(r') [ dv" dr" 



2-kTi 1 



V(r)V(r') / dr" 



V(r") 



V (r")V(r'") 



(B.20) 



from which we see that D(r',r) = D(r,r'). We will make use of these expressions in what 
follows. 

At present we are interested in the long- wavelength (low energy) properties of this reference 
T-matrix, hence we begin by expanding l(B.19() onto momentum eigenstates, yielding 



(klTfr, oik') = 



9 



dr dr' i 



-ik r ik'-r' 



H',V^I- ( 27r )3 

By applying the gradient operators Vk and V k we see that 



e 4k r £>(r,r'). 



(B.21) 



V k (k|T^|k') 



k,k'=0 



(2tt) 3 
-V*<k|2W|k' 



drdr' rD(r,r') 



lk,k'=0' 



(B.22) 



^From Eq. IB. 20(1 we see that J dr'D(r, r') is an even function of r provided only that V(— r) 
V(r). Hence the r.h.s of Eq. IB. 22(1 is the integral of an odd function which leads to 



v k (k|4,^|k')| kk , =0 = vwk|4, )t ,|k')| k!k)=0 = 0. 



(B.23) 



The fact that the gradient vanishes at k, k' = implies that this point is either an extremum 
or inflection point, i.e. the reference T-matrix is topologically flat in the long wavelength limit. 
Extending this approach will lead to the result that the curvature in k-space is proportional to 
the second moment of D(r, r') so that we can approximate the T-matrix by its long-wavelength 
limit provided that \k\R » |k'|i? -C 1, where R is roughly the radius of the kernel D(r,r'). 
Replacing the reference T-matrix by its long-wavelength limit and recalling that the kernel is 
normalized to unity gives 

(B.24) 



(k|IW|k'> 



9 

(27T) 



which then leads to the delta-function approximation 



<r|7W|r'} 



9 



(2tt) 6 
= 9S(r)5(v') 



(B.25) 



which is the main result of this section. 
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